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We report the first experimental demonstration of two- 
plioton imaging witli a pseudo-tliermal source. Similarly to 
the case of entangled states, a two-photon Gaussian thin lens 
equation is observed, indicating EPR type correlation in posi- 
tion. We introduce the concepts of two-photon coherent and 
two-photon incoherent imaging. The differences between the 
entangled and the thermal cases are explained in terms of 
these concepts. 

Two-photon imaging has been the subject of massive 
studies since its first demonstration ten years ago [1,2]. 
In a general sense two- photon imaging involves a joint de- 
tection of two photons at distant space-time points. The 
radiation from a source is split into two separate optical 
paths; an object (aperture) is placed in one of the paths, 
but the spatial information of the object is recuperated 
in a nonlocal fashion by means of the second order cor- 
relation measurement. The effect of two-photon imaging 
has been brought to the general attention by an exper- 
iment that exploited the quantum entanglement nature 
of a photon pair generated via Spontaneous Parametric 
Down Conversion (SPDC) [2]. In that experiment the 
signal and idler radiation of SPDC were sent to two dis- 
tant detectors. An aperture and an imaging lens were 
placed in the signal arm of the optical setup just before a 
bucket detector; there was no optical element in the idler 
arm; however, by scanning the idler detector in the plane 
defined by a two-photon Gaussian thin lens equation, a 
sharp and magnified image of the aperture was obtained 
in the coincidence counts, even though the single count- 
ing rates of both detectors were fairly constant during the 
scanning. This historical experiment was named ghost 
imaging due to the fact that the image of the object is 
formed by photons that never actually pass through the 
aperture. Further investigations on entangled photons 
brought to the development of a new field, named two- 
photon geometric optics [3]. 

Recently it has been argued that classically correlated 
light might mimic some features of the entangled pho- 
ton pairs in coincidence imaging setups [4]. Notice that 
the possibility of simulating the two-photon imaging fea- 
tures of entangled states with classical sources was not 
ruled out by the authors of the original ghost imaging ex- 
periment [2] . Both the theoretical work of Abourraddy 
et.al. [5] and the experimental investigation of Bennink 
et.al. [4] stimulated a very interesting debate about the 
role of entanglement in two-photon coincidence imaging 
[6-11]. In particular Bennink et al. [4] have experimen- 
tally demonstrated the possibility of performing "far- 
field" coincidence imaging, i.e., measuring momentum 
correlation only, with pairs of photons (pulses) classi- 
cally correlated in momentum. Very recently Gatti et.al. 
[12] have proposed thermal (or pseudo-thermal) radiation 



as a classical source to perform "near-field" coincidence 
imaging, which measures position correlation (at least 
partially), in a specific optical setup. Since then a great 
deal of attention on the subject has been induced [13]. 

In this paper we wish to present the first experimental 
demonstration of two-photon ghost imaging with thermal 
radiation. In particular, we show that a thermal source 
is able to simulate one of the main features of entangled 
two-photon imaging: the two-photon Gaussian thin lens 
equation, i.e., the EPR type correlation in position is 
partially observable. The expression partially accounts 
for the reduced visibility of thermal two-photon images 
(50% constant background). Furthermore, we introduce 
the concepts of two-photon incoherent and two-photon 
coherent imaging to explain the fundamental differences 
between thermal and entangled two-photon imaging. No- 
tice that the thermal ghost image presented here is the 
first authentic simulation of the quantum ghost image 
[2]. In fact, no two-photon thin lens equation, which 
is a symbol of position correlation, can be found by re- 
placing SPDC with a source of photons classically corre- 
lated in momentum [4] . Thermal sources provide an ideal 
comparison to quantum ghost imaging effects for a more 
complete understanding of two-photon, or two-particle 
physics. 

The behavior of entangled two-photon systems has 
been well studied [14,15]. It is possible to establish 
an analogy between classical optics and entangled two- 
photon optics: the two-photon probability amplitude 
plays in entangled two-photon processes the same role 
that the complex amplitude of the electric field plays in 
classical optics; the role played by the intensity of the 
electromagnetic field in classical optics is played by the 
rate of coincidence counts, and therefore by the time 
integrated second-order correlation function in entan- 
gled two-photon processes. For both thermal source and 
SPDC, the two-photon thin lens equation is the result of 
the coupling of two-photon probability amplitudes. How- 
ever, for an SPDC source, one pair of photons contains 
all the two-photon probability amplitudes that generate 
the ghost image; we define the result as two-photon co- 
herent image. For a thermal source, instead, the various 
two-photon probability amplitudes come from the ensem- 
ble of many independent pairs that generates the ghost 
image; thus, we define it as two-photon incoherent image. 

The experimental setup is shown in Fig. 1. After 
the pseudo-thermal source [16,17], a non-polarizing beam 
splitter {BS) splits the radiation in two distinct optical 
paths. In the reflected arm an object, with transmission 
function T{xi), is placed at a distance cIa from the BS 
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and a bucket detector (Di) is just behind the object. In 
the transmitted arm an imaging lens, with focal length 
/, is placed at a distance ds from the BS, and a multi- 
mode optical fiber (then connected to detector D2) scans 
the transverse plane at a distance d'^ from the lens. The 
output pulses from the two single photon counting detec- 
tors are then sent to an electronic coincidence circuit, to 
measure the rate of coincidence counts. 

The rate of coincidence counts is governed by the sec- 
ond order Glauber correlation function [18]: 

G(2)(ti,fi;f2,r2)^ (1) 

{e[-^ {h, n)Et^ {t2, f2)4+^ {t2, r2)E[-^^ (tun)). 

where E^~^ and £^(+) are the negative-frequency and 
the positive-frequency field operators describing the de- 
tection events at the space-time locations {f\,ti) and 
(i^, ^2)- The transverse second order correlation function 
(at equal times) for a thermal source is given by [19]: 

(^thermal (^1 ; ^2) OC ^ |fifi (f i , 1^ ^ IffS (^2 , |^ 

q q' 

+ I^<?^(^l,9)ff2(x2,9)|^ (2) 

where Xi is the transverse position of detector Di , q and 
q' are the transverse components of the momentum vec- 
tors and gi{xi;q) is the Green's function associated with 
the propagation of the field from the source to the i*^ 
detector [20]. Notice that there are two differences with 
respect to the SPDC case: (1) the presence of a back- 
ground noise (first term of Eq. 2), which does not exist 
for SPDC; (2) the possibility of writing the second term 
of Eq. 2 as a product of first order correlation functions, 
^12^^21^ while the second order correlation function for 
SPDC; cannot be reduced to a mere product of first order 
correlation functions in any format. 

For the setup of Fig. 1, it can be shown that for any 
values of the distances d^, ds, and d'^ which obey the 
equation: 



ds -dA d'g f 
Eq. 2 can be simplified as [19]: 

Gi^l(f2)cxiV+|T(^^f2)^ (4) 

where T^ C^'V'^^ ^2) is the object transmission function 
(T(xi)) reproduced on the D2 plane. We can then con- 
clude that a thermal source allows reproducing in co- 
incidence measurements the ghost image of an object, 
similarly to the SPDC case, except for a constant back- 
ground noise. The constant background noise, first term 
in Eq. 4, is proportional to the total transmittance of the 
object and, therefore, can be simplified to N, the number 



of transparent features in the object plane [19]. Eq. 3 can 
clearly be interpreted as a two-photon Gaussian thin lens 
equation, in which the object distance is So — ds — dA 
and the image distance is Sj = d'g (see Fig. 2). We ex- 
pect to observe an inverted image magnified by a factor 
of M = Si/so- Notice that when only one slit is inserted 
in the object plane, the maximum achievable visibility is 
33%. The visibility, however, is expected to drop when 
the number of transparent features increases. 

The discovery of general laws in physics allows con- 
sidering pictorial representations which may turn out to 
be powerful but still simple predictive tools. In this 
sense Klyshko's pictures for SPDC imaging experiments 
[21,2] are exemplar. The results presented in Eq. 3 and 
Eq. 4 also oH'er the possibility of considering a general- 
ized Klyshko's representation; nevertheless, the presence 
of two terms in Eq. 4 makes the thermal Klyshko's pic- 
ture more complex. For simplicity, in Fig. 2 we only 
picture the second term of Eq. 4, which is the one that 
produces the ghost image. Since the source is chaotic, 
each atom of the source randomly emits photons with all 
possible values of momentum. A photon with momentum 
q arriving at xi can produce a coincidence count with 
any other photon emitted by the source and arriving at 
X2- These coincidences result in a background noise that 
does not lead to the production of an image. This part 
of the physical process is not reflected in Fig. 2. Notice 
that, similarly to the SPDC case, the point-to-point cor- 
respondence between object and image plane is the result 
of the addition of two-photon probability amplitudes. It 
is important to point out however a remarkable differ- 
ence between the thermal and the SPDC case. In the 
Klyshko's picture for SPDC each pair of lines corresponds 
to a two-photon probability amplitude and all of them are 
associated to one pair, hence we define the image two- 
photon coherent. In the diagram for the thermal case 
each pair of lines represents the term gl{xi,q)g2{x2,q) 
and the various lines are associated to difi'erent pairs of 
photons. It is apparent then that the thin lens equation 
is the result of the coupling gl (xi , q)g2 {x2,q), but the 
image is generated by the ensemble of many independent 
pairs, leading us to the definition of two-photon incoher- 
ent image. The diagram of Fig. 2 also shows that while 
in the Klyshko's picture for SPDC the source behaves as 
an ordinary mirror, here the thermal source behaves as 
a phase conjugated mirror; this is due to the definition 
of the object distance Sg = ds — dA- This interpretation 
indicates the presence of a pseudo-object in the plane a, 
as shown in Fig. 2; scanning this plane the pseudo-object 
was actually observed (see also Fig. 3). 

Our first experimental measurement was aimed to ver- 
ify the existence of a point-to-point correspondence be- 
tween object and image plane, as expected by the ex- 
istence of a thin lens equation. The setup is the same 
as that of Fig. 1, but we used the 60/im-diameter input 
aperture of an optical fiber as the object (whose out- 
put was then coupled to detector Di). As a preliminary 
measurement we studied the temporal second order cor- 
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relation function. An optimized coincidence time window 
was chosen accordingly. For the actual spatial correlation 
measurement, we collected three sets of data for three dif- 
ferent positions of the fiber in the object plane; in every 
measurement we kept the position of the fiber fixed and 
scanned detector D2 in the transverse direction. The re- 
sults are shown in Fig. 3: any shift of the fiber in the 
object plane causes a shift of the correlation function in 
the opposite direction, in analogy to standard geomet- 
rical optics. In particular Fig. 3 shows that shifting the 
fiber in the object plane by 2mm the correlation function 
shifts by 4.3mm. Hence, the magnification in the imag- 
ing plane is Mmeas = 2.15, which is very close to the 
expected value {Mexpect — 2.16). The achieved visibil- 
ity is 26%. The results shown in Fig. 3 clearly show the 
point-to-point correspondence between object and imag- 
ing plane in agreement with Eq. 3 and Eq. 4. 

Then we placed a double slit in the object plane (cen- 
ter to center separation 1mm, slit width 0.2mm) and 
repeated the measurement. The results are shown in 
Fig. 4. As expected, the single counts are flat, while 
the coincidence counts reproduce the magnified image of 
the double slit. The visibility drops to about 12%. 

The experimental data show that the visibility of the 
two-photon image drops with the number of transparent 
features in the object plane, as predicted in Eq. 4. This 
effect is readily understood by inspecting Fig. 3: for each 
feature in the objc;c;t plane, the whole image plane shows 
besides the expected image a non-negligible noise level. 
Hence, if in the object plane there are simultaneously 
three features, in the image plane we will observe the ad- 
dition of the three graphs of Fig. 3, as predicted by Eq. 4. 
This clearly indicates that the background noise increases 
with the number of transparent features to image at the 
expenses of the visibility. 

The physics behind quantum ghost image and ther- 
mal ghost image can be understood as follows: in both 
cases the image is the result of the addition of two-photon 
probability amplitudes nevertheless such probability am- 
plitudes have very different origins in the two cases. In 
SPDC, a signal-idler pair is described by a pure state 
given by the superposition of an infinite number of two- 
photon probability amplitudes. On the other hand, ther- 
mal sources are modelled as an incoherent statistical mix- 
ture of many pairs of photons; the various two- photon 
probability amplitudes are provided by the entire ensem- 
ble of photon pairs. The different models come from the 
different processes involved in the generation of SPDC 
and thermal radiation. SPDC comes from a nonlinear 
coherent interaction and consequently it is impossible to 
identify the birth place of a photon pair. This is why 
each photon pair contains simultaneously all the possible 
transverse momenta necessary for imaging. On the other 
hand, for a thermal source each photon is generated from 
a random excitation of independent atoms. This is in 
line with the interpretation of the two Klyshko's pictures 
given above and clarifies what we mean by two-photon 
coherent and two-photon incoherent imaging. 



In conclusion, we have presented the first experimen- 
tal demonstration of two-photon ghost imaging with 
thermal-like sources. For the first time a two-photon 
Gaussian thin lens equation is found for non-entangled 
sources, indicating the existence of, at least partial, EPR 
type correlation in position. We have also introduced 
the concepts of two-photon coherent and two-photon in- 
coherent imaging to describe the different physics behind 
entangled and thermal ghost imaging. The already poor 
visibility of the thermal ghost image (max 33%), drops 
rapidly for complicated apertures; therefore, practical ap- 
plications of thermal ghost imaging rely on the develop- 
ment of a detection scheme which is insensitive to the 
background noise. Such detection scheme is now under 
development in our laboratory. 
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Figure 2. Alejandra Valencia, Giuliano Scarcelli, 
Milena D'Angelo, and Yanhua Shih. 
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FIG. 1. Experimental setup. Source diameter ~ 200/xm; 
a = 125mm; dA = 88mm; ds — 212mm; d'g = 268.5mm 
/ — 85mm; fiber tip diameter — 60^. 



FIG. 2. Conceptual unfolded version of the optical setup. 
Object plane and image plane obey a two-photon gaussian 
thin lens equation and are defined by So = ds — dA and 
Si = d'g, respectively. In terms of Klyshko's picture, the 
thermal source acts as a phase conjugated mirror forming a 
pseudo-object in the a plane. 



FIG. 3. Normalized second order correlation function vs. 
position of D2. The measurement shows the point-to-point 
correspondence between object and image planes, a) Tip of 
the fiber in the object plane located in the position indicated 
by the square, b) Tip of the fiber in the object plane located 
in the position indicated by the circle (central position), c) 
Tip of the fiber in the object plane located in the position 
indicated by the triangle. 



FIG. 4. Single and coincidence counts vs. transverse posi- 
tion of D2 in the image plane (2:2). Single counts of both Di 
(hollow circles) and D2 (full circles) are fiat. The coincidence 
counts (solid line with circles) show a magnified image of the 
object. 
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Figure 3. Alejandra Valencia, Giuliano Scarcelli, 
Milena D'Angelo, and Yanhua Shih. 
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Figure 4. Alejandra Valencia, Giuliano Scarcelli, 
Milena D'Angelo, and Yanhua Shih. 



Figure 1. Alejandra Valencia, Giuliano Scarcelli, 
Milena D'Angelo, and Yanhua Shih. 
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